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Abstract 



We report the first definitive experimental observation of periodic orbits 

(POs) in the spectral properties of an elasto dynamic system. The Fourier 

transform of the density of flexural modes show peaks that correspond to 

stable and unstable POs of a clover shaped quartz plate. We change the 

shape of the plate and find that the peaks corresponding to the POs that 

hit only the unperturbed sides are unchanged proving the correspondence. 

However, an exact match to the length of the main POs could be made only 

after a small rescaling of the experimental results. Statistical analysis of the 

level dynamics also shows the effect of the stable POs. 
PACS number(s): 05.45. +b,05.45.Mt,46.40.-f 
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Mechanical wave systems show locahzation and band gaps similar to solid state systems 
raising the possibility of applications An issue of current interest is if techniques 

developed in the analysis of quantum wave systems with chaotic classical counterparts can 
be applied to elastodynamic systems. The first experiments on elastodynamic systems were 
conducted with aluminum blocks and a statistical analysis was carried out . Waves can 
propagate in transverse and longitudinal modes in a solid and the application of analysis 
techniques developed in the context of quantum systems is not straightforward. There are 
two main approaches to analyzing quantum wave systems. The first is a statistical approach 
where the properties of the eigenvalues can be compared to the universal predictions of 
random matrix theory (RMT) PJ^. The second semiclassical approach uses the periodic 
orbits (POs) in the system as a framework to understand the structure of the eigenvalues [0 . 
In recent experiments it has been shown that the flexural modes of plates show statistical 
properties described by RMT [P-[TO||. However, the experimental observation of POs in 
elastodynamic systems has proven to be difficult. 

A semiclassical theory using the biharmonic equation for the flexural modes of a thin 
plate suggests that a close correspondence between the modes and ray orbits may be expected 
for the elastodynamic system because of the similarity of the governing wave equations []TT] . 



However, experiments indicate that higher order corrections are required to capture the 
dispersion relation in plates at higher frequencies ||10|- Therefore the biharmonic equation 
is not always a good approximation for the flexural modes of a plate in a regime where 
semiclassical analysis can be applied. We reported in an earlier study using a clover shaped 
plate that the statistical properties indeed show deviations from universality due to the 
presence of stable POs, but the result of the periodic orbit analysis was unclear Earlier 



studies with aluminum blocks have also yielded ambiguous results for the periodic orbits [jT3 

In this paper, we report the flrst detailed experimental study of the signatures of the 
periodic orbits in the spectral properties of an elastodjTiamical system. We perform mea- 
surements on a clover shaped quartz plate which has both stable and unstable POs (see 
Table I). Fourier transform is performed on the flexural density of modes using recently 



developed theoretical results for the dispersion relation [|T0|- To test the robustness of the 



peaks, we change one of the parameters which defines the shape of the plate in small steps 
that also alters some of the POs. The peaks corresponding to POs which do not hit the 
perturbed side are observed to remain unchanged. Furthermore, the peaks corresponding to 
the POs which hit the perturbed side are observed to change, proving the correspondence. 
An exact match to the length of the main POs could be made only after a small rescaling of 
the experimental results. We also analyze the statistical properties of the elastodynamical 
resonances under parametric variation using the parametric number variance v{x) and the 
velocity autocorrelation c(x). 

In the experiments, a fused quartz plate of thickness h = 1.5875 mm with dimensions 
given in the caption to Table I is used. The plate is kept on three piezoelectric transducers 
(one transmitter and two receivers) and the vibrations of the plate are measured using a 
HP4395A Network Analyzer. Therefore the plate vibrates freely inside the chamber, the 
only contacts are through tiny ruby spheres that are attached to the transducers which do 
not perturb the resonances. The plate and the transducers are placed inside a temperature 
controlled chamber at 300 K and at a pressure below 10^^ Torr to prevent losses due to air 
damping. A more detailed description of the experimental setup can be found in Ref. 0. 

The vibrational modes of a thin plate can be divided into fiexural and in-plane modes. 
These two types of modes obey different dispersion relations and it is therefore important 
to separate them before analyzing the data. To experimentally isolate the fiexural class of 
modes, we increase the pressure inside the chamber, and measure the resonance spectrum at 
~ 300 Torr. The Q-factor of the fiexural modes which is about 10^ decreases, whereas the 
Q-factor for the in-plane modes are unchanged P, p!2| , p!^ . Thus the resonances corresponding 
to the fiexural modes are identified. 

In the first series of experiments, we sanded off material at one of the edges of the plate 
[indicated by thick line in the inset to Fig. |l](a)] in 61 small steps. This causes PO number 
3 (shown in Table I) to decrease in length while PO number 4 remains unchanged. The 
fiexural modes were measured for each step in a frequency interval between 52 kHz (mode 



number 80) and 352 kHz (mode number 605) with a frequency resolution of 5/8 Hz. We 
note that all modes in this frequency regime are two dimensional as the wavelength is much 
greater than the thickness of the plate. In the second series of experiments, the same edge 
is first sanded at an angle to completely destroy PO number 3 and the resulting shape is 
shown in the inset to Fig. |I](b). Then the same edge is sanded in 60 small steps and the 
flexural modes in the frequency interval between 52 kHz (mode number 78) and 252 kHz 
(mode number 398) are measured. 

The change in the flexural modes as a function of mass removed m due to sanding is 
shown in Fig. |I|for the two series of experiments (the total mass of the unperturbed plate was 
21.35140 g). We observe that a number of modes are observed to split in Fig. |l](a) because 
the C^y point group symmetry of the shape of the plate is broken. One may expect some 
modes that are quantized about unperturbed PO number 4 to remain unchanged. However, 
none of the modes are observed to remain constant. 

According to theoretical expectations, |-F(OP versus length / should show strong peaks at 
values corresponding to the length of the POs of the ray system which are listed in Table |. 
Therefore to check for the signatures of the POs in the elastodynamical spectra, we calculate 
the square of the Gaussian-weighted Fourier transform for all the 121 sets of flexural modes 
using 
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(1) 



where p{k) represents the density of modes J2i<n^{k — hi), n is the total number of modes 
in a spectrum, and ki denotes the wavenumber of the iih mode. The parameter j3 is chosen 
such as e~^'^™'^== = 0.01, a value sufficiently small to reduce the spurious oscillations of the 



peaks due to finite size effects, and r = following Ref. [11 



A new dispersion relation that relates the wavenumber k to the frequency / for plates 



has been calculated in Ref. |10] and is given by: 



12^/^ n- 

k = —^VVL{l + aiVL + a2^^), (2) 



where n = y2/(l — iy), f2 = 2nfh/cs, and u is the Poisson's ratio. For fused quartz, 
u = 0.16, Cs = 3750 ms~^, ai = 0.177, and a2 = —9.40 x 10^^. Until these recent results, 
only the first term in the expansion given in Eq. ^ was used to calculate k. 

The P was obtained for each of the 121 sets of data and is plotted against 0.92 x / in 
Fig. ^. A number of strong peaks are observed. The locations and the strength of some of the 
peaks change while others remain constant. As noted earlier, the first series of perturbation 
decreases the length of the PO number 3 while keeping the PO number 4 constant. Indeed 
we observe in Fig. ||(a) that the peak at the length corresponding to these orbits splits with 
one peak shifting to a lower value and the other remaining unchanged (see Inset to Fig. |^ 
for detail). After the plate is perturbed to destroy PO number 3, the corresponding peak in 
Fig. 0(b) is observed to disappear. 

We note that even the amplitude of the peaks changes in a way which is in agreement 
with the idea that the area of the tori associated with that PO contributes to the amplitude 
of the peak. For example the amplitude of the peak at 20.3 cm is decreasing as the plate is 
sanded and one of the tori is slowly destroyed (see Fig. |^). The amplitude is constant in the 
second series of experiments after the PO number 3 is completely destroyed (see Fig. ||(b)). 

However while making a comparison with the actual length of the POs, we find that the 
peaks in the |-F(OP occur at a slightly higher value. This is the reason for the factor 0.92 
which multiplies the length I while presenting the results in Fig. The source of this factor 
is unclear and does not appear to correspond to errors in the values of u and Cg. After the 
scaling, several other peaks can be matched with the POs listed in Table |. Furthermore, 
it can be noted that the peaks corresponding to the POs which hit only the unperturbed 
dimensions are observed to be unchanged, while others are observed to change. 

We note that because the peaks are somewhat shifted our earlier study in which we 
examined only the |-F(/)P for one set of data was inconclusive. By parametrically varying 
the modes and examining the Fourier transform of each set of modes, we are able to now say 
with confidence that the peaks in the |-F(OP correspond to the POs of the elastodynamic 
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system. 

We also note that if the k is calculated using only the first term in the dispersion relation 
given in Eq. 0, then the correspondence is not observed. Therefore the corrections to the 
dispersion relation which takes into account the three dimensional nature of the plate are 
necessary to make the correspondence. It is possible that an addition correction in the 
theoretical formula is required to account for the slight scaling which is still required to get 
an exact match. 

We further analyze the statistical properties of the eigenvalues under parametric variation 
using the parametric number variance v{x) and the velocity autocorrelation c{x). Deviations 
may be expected from the universal forms for chaotic systems calculated using the 0-d 
nonlinear cr-model [jl5[ because of the presence of stable POs. Most systems in nature have 



both stable and unstable POs, and therefore it is interesting to know the nature of the 
deviations. Therefore we also investigate the properties of the eigenvalues under parametric 
variation. The phase space explored by ray trajectories inside the clover geometry has been 
discussed in detail in Ref . . The clover geometry has both chaotic and integrable regions 
in its phase space and it is therefore an example of a mixed system. To analyze the data, 
we normalize the energy E corresponding to the flexural modes in units of the local mean 
level spacing A to obtain the normalized energy e = E/A and the parameter m in units 
of the square root of the local mean squared slope to give the rescaled external parameter 

The parametric number variance v{x) which measures the collective motion of levels 
under parametric variation [jl^ is defined as v{x) = {{n{e,x') — n{e,x' + x))^) , where 
n{e, x) is the staircase function which counts the number of energy levels at fixed x with 
energy lower than e, the average is over the normalized parameter x' and e. v{x) for the 
clover shaped plate with C^^j point group symmetry and the plate without any symmetries 
is plotted in Fig. §(a). v{x) for the plate without symmetries is observed to be greater than 
the universal form for a completely chaotic system |T3|. Greater deviations are present in 



case of the data corresponding the first series of experiments, where the symmetry of the 
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plate is being broken and degeneracies lifted. 

To study the oscillations of the eigenvalues during perturbation, we calculated the ve- 
locity autocorrelation c{x) = where the average is over the parameter x' 



dx' dx' 

and the energy e. For this correlator no analytical results exist for the intermediate values 



of X therefore we compare our results to a curve numerically calculated by Mucciolo |[T8 
which agrees with the analytical results in the limit of large and small x |]15|. This form 
also agrees with the correlations observed in the spectra of plates with shapes that are com- 
pletely chaotic The result for c{x) for the two series of experiments is shown in Fig. 0(b). 
The curve corresponding to the clover shaped plate without symmetries deviates from the 
universal results similar to that of v{x). Therefore the presence of stable POs makes the 
oscillations slightly more irregular compared to the completely chaotic case which leads to 
smaller anti-correlations. Because of the added effect of breaking symmetry, the deviations 
from the universal results are even greater in the first series corresponding to Fig. |l|(a). 

In summary, the study provides the first conclusive evidence of the relevance of the 
periodic orbits in the analysis of elastodynamic systems. The ability to make this corre- 
spondence was possible only because of the systematic changes in the peaks of the Fourier 
transform due to shape deformation. However, even after using the improved relation for 
the dispersion relation, small discrepancy between the actual location of the peaks and the 
length of the periodic orbits is observed. We have further shown that the presence of the 
stable periodic orbits causes the level motions to fluctuate differently than the universal 
forms for completely chaotic systems. 

We thank Olivier Brodier, Kristian Schaadt, and Taimur Ellahi for assistance in setting 
up the experiments, and Joel Norton for technical assistance. This work was supported by 
Research Corporation. A.K. was partially supported by a Alfred P. Sloan Fellowship. 
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TABLES 

TABLE I. The periodic orbit (PO) and their length Ipg for the clover shaped plate. The radius 
of the convex side is 3.556 cm, the radius of the concave side is 5.080 cm, and the distance between 
the longitudinal concave sides is 10.160 cm. PO number 1, 2, 5, 6 and 8 are unstable and the rest 
are stable. 
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FIGURES 

FIG. 1. (a) Representative motion of the unfolded flexural modes of the clover shaped plate 
as a function of the removed mass m obtained by sanding the edge of the plate indicated by the 
thicker line in the inset. The perturbation preserves the reflection symmetry about the vertical 
axis and decreases the length and stability of PO number 3. (b) The motion of the unfolded level 
after the reflection symmetry of the plate and PO number 3 is destroyed. 

FIG. 2. (a) versus length / rescaled by a factor of 0.92 for the first experimental series. 

|F(Z)|2 is obtained using approximately 500 flexural modes in the frequency interval 52 kHz to 352 
kHz. (b) versus length I rescaled by a factor of 0.92 for the second series, is obtained 

in this case using 306 flexural modes in the frequency interval 52 kHz to 252 kHz. The position 
of the POs of the system are also indicated. The numbers corresponds to the POs shown in the 
Table. |. 

FIG. 3. The ratio of the amplitude A of the peak near PO number 4 normalized by the ampli- 
tude of peak Ai corresponding to the 1st data set. Inset: A detailed plot of the peaks in |-F(OP of 
the 1st, 31st and 61st data set in the parameter m that correspond to PO number 3 and 4. 

FIG. 4. (a) The parametric number variance v{x) versus the normalized parametric variation 
X due to sanded mass m. (b) The velocity autocorrelation c{x) versus x for the experimental data 
compared to universal RMT curves. The deviations in the case of the second series is due to the 
presence of stable POs. The deviations are stronger in the first data series because the symmetry 
of the plate is being reduced. 
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This figure "Figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0110543vl 



TABLES 



Number PO /, 



6 



po 

(cm) 









14.5 



16.7 



20.3 



20.3 



20.5 



21.1 



Number PO /, 



10 



11 



12 



po 

(cm) 









22.4 



23.4 



28.7 



31.2 



38.9 



64.5 



TABLE I. The periodic orbit (PO) and their length Ipo for the clover shaped plate. The radius 
of the convex side is 3.556 cm, the radius of the concave side is 5.080 cm, and the distance between 
the longitudinal concave sides is 10.160 cm. PO number 1, 2, 5, 6 and 8 are unstable and the rest 
are stable. 
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This figure "Fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0110543vl 
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FIG. 3. The ratio of the ampHtude A of the peak near PO number 4 normahzed by the amph- 
tude of peak Ai corresponding to the 1st data set. Inset: A detailed plot of the peaks in of 
the 1st, 31st and 61st data set in the parameter m that correspond to PO number 3 and 4. 
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FIG. 4. (a) The parametric number variance v{x) versus the normahzed parametric variation 
X due to sanded mass m. (b) The velocity autocorrelation c(x) versus x for the experimental data 
compared to universal RMT curves. The deviations in the case of the second scries is due to the 
presence of stable POs. The deviations are stronger in the first data series because the symmetry 
of the plate is being reduced. 
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